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Contact electrification of identical insulating particles has crucial significance for industrial and 
environmental science, especially in wind-blown granular systems. At the same time, the experimental 
phenomena of charge transfer first increased and then decreased with the increase of relative humidity has 
attracted the interest of many researchers. Humidity can affect the charge transfer has been early observed in 
the experiment, but the reason always puzzles researchers. In this study, based on trapped high-energy 
electron transfer mechanism, we introduce the effect of the water film in the charge transfer model and 
consider the actual situations of the sand particles in the collision process. Furthermore, charge transfer 
between sand particles in a single collision under different humidity conditions is investigated. The 
predicted results agree well with the law obtained in existing experiments qualitatively and thereby a 
possible explanation why humidity can affect the charge transfer is given. 



The electrification of sand is a natural phenomenon in wind-blown granular systems, e.g., dust devil, sand/ 
dust storm and wind-blown saltation. The earliest observation on electrification of sand can be traced back 
to 1914 conducted by British physicist Rudge 1 who performed a measurement experiment of the electric 
field in the sand dust storm. This contact electrification of identical materials phenomenon also exists in other 
granular systems, including fluidized beds 2 , powder handling 3,4 , volcanic eruptions 5,6 , wind-blown snow 7 . In 
wind-blown granular systems, both the measurements and the theoretical predictions of the polarity of charged 
sand particles indicate that smaller particles generally charge negatively, while larger particles charge posi- 
tively 8 " 12 . Although existing theoretical models can explain the contact electrification of identical materials of 
part phenomenon observed in experiments to some extent, it's still not perfect. And we know little about the 
microscopic mechanism of how humidity can affect the charge transfer. 

Humidity can affect the charge transfer has been observed in many experiments. The electric field changes with 
the relative humidity had been early observed by Lowell and Truscott 13 and Pence et al. 14 . Guardiola et al. 15 
adopted a fluidized bed of glass beads experiment to investigate the influence of humidity on the degree of 
electrification and indicated that the variance of the charge with humidity existing five regions, and first increase 
then decrease to zero with the humidity. The wind tunnel experiment of Qu et al. 16 indicated that moisture 
content could influent the electric field. In the experiment, the electric field intensity increases with moisture 
content in the region of 0.5% ~ 1.0%, but when it increases to 2%, electric field intensity decreases, which suggests 
that humidity can influent the charge transfer. Ducati et al. 17 found that the change of environment humidity can 
make the insulator surface charge easily. The wind tunnel experiment of Xie and Han 18 showed that aeolian 
electrical field intensity linearly increases when the relative humidity is less than critical relative humidity, but 
exponentially decreases when the relative humidity is higher than the critical value. And the critical relative 
humidity is affected by the wind speed. Giffm and Mehrani 19 measured the net charge of three regions of the 
fluidized bed at two fluidization gas velocities and different relative humidities, and they directed that humidity 
can affect the charge in all regions and the variations of charge with humidity are not unique. 

Although the above experimental results have proven that water plays a significant role in contact electrifica- 
tion process of identical insulators, how humidity can affect the charge transfer so far still has not obtained a good 
explanation within the theoretical simulation. The insulator contact electrification mechanism mainly has two 
viewpoints, one is based on the mechanism of ion transfer, and another is based on the mechanism of electrons or 
holes transfer. 

The representative work for the ion transfer is completed by Latham et al. 20 . They regarded that there exists a 
temperature contrast, between the sand particles during the collision due to the difference of the sand surface area. 
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This temperature contrast leads to the certain ions (H+ or OH — ) 
transfer along the direction of the temperature gradient and then 
sand particles charge. Subsequently, other scholars 17,21 " 23 carried 
out the further researches based on the ion transfer. However, it is 
difficult to measure the temperature contrast and the charge between 
two sand particles of different diameters during the collision. 
Therefore, it is hard to be proved from the experiment. And 
researches have shown that the identical insulator surface can be 
charged even in the condition of absence of water 13 . Thus, the ion 
transfer model cannot describe the contact electrification phenom- 
enon in a dry environment. 

Lowell and Truscott 24 suggested the asymmetric contact electri- 
fication theory of identical insulators based on high energy trapped 
electrons. And many scholars 25 " 29 expanded and optimized this the- 
ory to further explain the causes of contact electrification and pro- 
cesses of sand particles or other insulating particles in the last five 
years. Lacks and Levandosky 26 proposed a simple charge model to 
explain the tendency of smaller particles charge negatively while 
larger particles charge positively based on the explanation of 
Lowell and Truskott 24 . Kok 28 further simplified the trapped high- 
energy electrons transfer rule and assumed that electrons on the 
surface of particles can be in either a low-energy or high-energy state 
and electrons can be only released from high- energy states to low- 
energy states. Besides, he assumed that the amount of high-energy 
electrons transfer depend on particle diameters because of the exist- 
ence of maximum distance that an electron can tunnel during a 
collision. Hu 30 considered different contact areas using the elastic- 
plastic contact theory instead of the point contact in Kok's model 28 
and thought that holes release lead to the insulator charging. 
However, electrons or holes transfer mechanism also cannot explain 
how humidity can affect the charge transfer. 

Generally speaking, the mechanism of ion transfer is based on the 
temperature contrast which is hard to be measured between two 
different diameter sand particles during the collision. And although 
the model of Hu 30 is more close to the actual situation which con- 
sidered the plastic deformation, when adding the water film, using 
this model will increase the difficulty and complexity, besides, the 
cracking of water film is not considered, which means the surface of 
particles do not contacted in the collision process. Therefore, in this 
paper, the electron transfer model of Kok 28 is adopted to discuss the 
effect of the water film. The nanoscale pores of the sand surfaces 31 
and extrusion deformation of the water film are considered in our 
model. 

We introduce the physical model which considers the thickness of 
the water film, see Fig. 1. Then we can get the radius of the sand 
particle i and; which covered the water film: 



R , i = R i + n i , R , j = R j + n j . 



(1) 



Kok 28 assumed that electrons on the surface of particles can be in 
either a low-energy or a high-energy state. The numbers of high and 
low energy electrons on particle i and j at time t are represented as 
«/h(£)> n iL (t), rij H (t), rij L (f). The density of high-energy surface states 
are both represented as p 32 33 . Then the initial numbers of high- 
energy states can be obtained: 



n iH (0) = 4nR?p, n jH (0) = 4nR j 2 p. 



(2) 



Kok 28 also assumed that electrons can only release from high-energy 
states to low-energy states and cannot release to its own. For particle 
/, the net charge after one collision at time t is 

A qi (t) = e[n iH (0) - n iH (t) + n iL (0) - n iL (t)} 



(e is theelementary charge). 



(3) 




Figure 1 | Schematic of the charge transfer occurring during a collision 
between two different diameter sand particles which are covered by the 
water film. R{ and Rj are the radius of sand particle i and j y tii and rij are the 
water film thickness of sand particle i and j (rij = rij in the simulation), and 
3 0 (n) is the maximum distance that an electron can tunnel between two 
sand particles. 

with the numbers of low-energy states of particle j obtain from 
high-energy states of particle i, i.e., 

- An iH = An jL = InpR^l - cos tp t ). (4) 

In our model, we can get the following algebraic relations: 

Rf sin (p i = (Rj + d 0 (n))smoCj, (5) 



R'i + R) = (Rj + ^(w^cosa,- + R,- cos <p t . (6) 
From equation (5) and (6), we can derive the following relation: 
<5 0 2 («) + 2R;<5 0 («) - IRjto + nj) - (n^n-f 



1 — cos (p- = , , 

Tl 2R f (R' f + R' ; -) 

Then we can rewrite equation (4): 



-• (?) 



5 0 2 (n) + 2Rj5 0 (n) - 2R,-(n,- + «,■) - (», + «,) 2 
(R'i + R'i) 

,V(n) + 2R,A,(n) - 2R,-(n / + n,-) - (n,- + n,-) 2 
(R'; + R'i) 

Finally, the relationship of charge transfer can be described as: 

(( 5 0 2 (m) _ (m . + m . ) 2 )(R ._ R . ) 



— AriiH = Ariji = npRj 
As symmetry, 

— AnjH=Ann = npRj 



Aqi = e[— AniH — Ann] = xep - 



(R'/ + R';) 



(8) 



(9) 



:io) 



And we can know that the numbers of high-energy state electrons 
of particle i release to low- energy states of particle j are accordant 



Compared with the results of Kok 28 , our theoretical equation (10) 
cannot only reflect the tendency of smaller particles charge nega- 
tively while larger particles charge positively but can also take the 
effect of the thickness of water film into account. Because n t and rij are 
too small comparing with grain radius Ri and Rj, so their effects on 
the R'i and R'j can be neglected. We all know that the thickness of a 
water molecule is about 0.4 nm. And the maximum distance S 0 (n) is 
also in the nanometer magnitude 28 . Then we will analyze the max- 
imum distance S 0 (n) quantitatively. 

The existing literature shows that, at the temperature of T = 23 °C, 
the surface conductivity a increases with the increase of the water 
film layers absorbed on the surface of quartz 34 (the most components 
of sand particles), and satisfies the following relations: 
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Table 1 | th 


3 relationship 


between 


relative humidity 


and absorbed water film 


layers 












Water film 
layers 


0 1 


2 


3 4 


5 6 


7 


8 


9 


10 


11 


12 


Relative 
humidity(%) 


0 4 


8.5 


20.9 31.5 


37.3 43.7 


48.5 


53 


56.9 


60.1 


64.1 


68.1 



3.0xl0- 18 xl0°- 44n («>1) 
6.5 xlO- 18 O = 0) 



(11) 



Where n is the absorbed water film layer. The relationship between 
surface conductivity a and energy barrier E b described as follows 35 : 



E b (n) = -2KT]n( — ) 



(12) 



Where k is the Boltzmann constant, T is thermodynamic temperature. 
We can see that the increase of the surface conductivity makes elec- 
trons through the energy barrier easier. The relation between energy 
barrier E b and the maximum distance S 0 (n) can be described as 36 



<5o(») = 



h 



sJ%mE b {n) nm a 



In ( — — ) 



(13) 



Where h is the reduced Planck constant, a is the radius of the well 
corresponding to the electron trap, t co u is the time scale of the collision, 

1.12 



m is the electron mass, and rj = 



- (2 + 71). Thus we can obtain the 



relational expression between the maximum distance d 0 (n) and the 
absorbed water film layers: 



So(n) = 



h 



\6mKT\n 



ln( — (14) 



In the equation (14), the maximum distance S 0 (n) increased with the 
increase of the absorbed water film layers. Here, we notice that the 
item (d 0 2 (n) — (tii + nj) 2 ) in the equation (10) will become interesting 
and may get the results that we want. 



Results 

It is worth mentioning that we also consider the actual situations of 
the sand particles in the collision process, (z) The nanoscale pores 
exist in the sand particle surface 31 and can absorb water molecules in 
the air. Based on this, when particle begin to absorb some layers of 
water molecules, water molecules do not immediately cover the sand 
particle surface to form the water film but exist in the nanoscale 
pores. It does not increase the collision distance but increase the 
conductivity until the water molecules fill the nanoscale pores. 
After that, the water film will be formed on the sand particle surface. 
And then this part («,■ + nj) 2 will be calculated by the equation (10). 
(zz) The thickness of the water film is not invariant during the col- 
lision. Because of the extrusion during the collisions, the thickness of 
water film will be compressed. Here, we define a compression ratio as 



compression value 
initial thickness of water film 



(15) 



Summer and Pitts 37 gave the relationship between absorbed water 
film layers on the glass surface and the relative humidity. Their 
results showed that the absorbed water film layers increase with 
the increase of the relative humidity. It is very difficult to directly 
measure the thickness of the water film on a tiny sand particle sur- 
face. So in our simulation, we use the data of glass which the com- 
positions close to the sand. In order to be more intuitive and easy to 
compare, we use the relative humidity to replace the absorbed water 
film layers in our results, and the corresponding relations shown in 
table 1. 

First, we consider the existence of the nanoscale pores on the 
surface 31 to investigate the effect of the thickness of water film on 
charge transfer. And we study whether the water film can be formed 
when the sand particle absorbs N ab =1 — 5 layers of water molecules. 
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Here, we adopt the energy barrier is 4.5 ev 28 when there is no water, 
for the given physical model, only the item of (3o 2 (n) — (n\ + nj) 2 ) 
changes in the equation (10), so we discuss the change of this item to 
reflect the influence of humidity on the charge transfer. And when 
the value of (d 0 2 (n) — (ni + nj) 2 ) is negative, it means an electron is 
no longer able to tunnel during the collision rather than transfer 
reversely. 

Fig. 2 shows the variation of this item with humidity at different 
absorbed water layers in a single collision. We can see that when the 
more absorbed water molecules needed to form the water film, the 
larger the critical humidity of charge transfers. As we can see, in N ab 
= 5 case, the critical humidity of charge transfer is about 43.7%, and 
the charge transfer first increases then decreases to zero, which is 
accordant with the experimental results of Guardiola et al. 15 . In spite 
that the measurement results of Guardiola et al. 15 are macroscopical 
quantities and particles average, while our results just show the 



charge transfer between two sand particles, but our results also can 
reflect the rules from a certain extent. In addition, in the results of 
Guardiola et al. 15 , we can see the influence of the relative humidity on 
the charge is considerable, especially the increment at the peak is 
obvious (the minimum increment of the given data is 21%), but in 
this model, it is small. 

Hence, in order to keep accordant with the actual situation, we 
improve our model by considering the compression of the water film. 
Reference on the former step, we assume that the sand particle 
becomes to form the film when it absorbs 4 layers of water molecules. 
Fig. 3, Fig. 4 and Fig. 5 shows the variation of charge transfer with 
relative humidity in a single collision when the energy barrier are 4.5, 
2.5, 1.0 ev respectively (The different components of the sand part- 
icles may result in the energy barrier variation). We can see that in the 
three cases, along with the increase of compression, the charge trans- 
fer experience three situations (a)first increase then decrease sharply, 
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Figure 4 | Charge transfer varies with relative humidity at different compressions when the energy barrier is 2.5 ev. 
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Relative Humidity(%) 

Figure 5 | Charge transfer varies with relative humidity at different compressions when the energy barrier is 1.0 ev. 



(b) first increase then decrease slowly, (c)only increase. From Fig. 3, 
when analyze the data at the compression of 88%, we find that the 
peak of the charge is at 35% humidity, and when the humidity is at 
70%, the charge transfer tends towards zero, which agree with the 
results of Guardiola et al. 15 . Compare Fig. 3, Fig. 4 and Fig. 5, we can 
find that reducing the energy barrier can effectively increase the 
charge transfer, and make the peak increments increase significantly. 
From Fig. 5, when analysis the data at the compression of 60%, we 
find that the peak of the charge is at 40% humidity, and the peak 
increment is 26% (compare to the humidity is zero). Besides, we find 
that the increment is much more obvious after considering the com- 
pression, and the tendencies of most curves are accordant with 
Guardiola et al. 15 , only a fraction of curves consistently increase with 
relative humidity. While Giffin and Mehrani 19 indicated that the 
variations of net charge with humidity obey different rules for dif- 
ferent measurement locations in their fluidized bed experiment. One 
of his data attracts our attention, which showed that at a large flu- 
idization gas velocity, the charge constantly increases with humidity, 
and it is accordant with part of our results. Our results are derived at 
large compression, and Giffin 19 is derived at large fluidization gas 
velocity, so the coincidence can be explained as: when the velocity of 
the sand particle is large, the contact is more violent, and compres- 
sion of the water film is much larger. The law of the experiment of 
Giffin 19 is accordance with Guardiola et al. 15 at small velocity, which 
is embodied in our simulation. At the same time, we note that the 
corresponding relative humidity at the peak of charge increases with 
the fluidized bed velocity of Guardiola et al. 15 , which is the same with 
our results. 

Discussion 

In summary, based on trapped high-energy electrons transfer theory, 
considering the actual situations of the sand particles in collision 
process, a possible theoretical model is proposed. We establish the 
relationship between relative humidity and charge transfer and 
simulate the effects qualitatively. The results show that our simu- 
lation results are accordant with the universal experimental results of 
the charge transfer first increased and then decreased with the 
increase in relative humidity. Our analyses also suggest that the 
surface conductivity increases with the relative humidity increasing. 
Thus, the trapped high-energy electrons can tunnel easier at first, 
then harder and decrease to zero after the relative humidity reaches a 



critical humidity. In particular, when the water film is not cracking, 
we can speculate that if the velocity of sand grain increases, the 
compression of the water film will become much severer accordingly. 
However, the water film may rupture at a critical velocity, and then 
the charge transfer decreases with this part of surface conductivity 
decreasing. It may be a reason why the effect of wind velocity is the 
same as humidity on the aeolian electric fields. 

In actual industrial production, the contact electrification of ident- 
ical insulating particles is a common problem. Through above 
researches, we can control not only the humidity but also the dimen- 
sion, surface condition and velocity to strengthen or weaken the 
charge transfer. 

Methods 

We establish a charge transfer physical model when the sand particle surface is 
covered with the water films. Then we consider the actual situations of the sand 
particles in the collision process and parameterize the assumptions. Finally, we carry 
out a numerical simulation of different hypothetical cases and compare with the 
existing experimental results. 
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